Unified description of the dynamics of quintessential scalar fields 
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Using the dynamical system approach, we describe the general dynamics of cosmological scalar 
fields in terms of critical points and heteroclinic lines. It is found that critical points describe the 
initial and final states of the scalar field dynamics, but that heteroclinic lines give a more complete 
description of the evolution in between the critical points. In particular, the heteroclinic line that 
departs from the (saddle) critical point of perfect fluid-domination is the representative path in 
phase space of quintessence fields that may be viable dark energy candidates. We also discuss the 
attractor properties of the heteroclinic lines, and their importance for the description of thawing 
and freezing fields. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The dynamics of scalar fields has been under a careful 
scrutiny since their involvement as strong candidates in 
dark energy models; those models useful for dark energy 
are generically dubbed as quintessence fields, in order to 
distinguish them from purely inflationary scalar fields 0- 
3]. There seems to be a consensus in that there are two 
types of extreme behaviors for quintessence fields: thaw- 
ing and freezing Thawing fields have at the beginning 
a behavior quite similar to that of a cosmological con- 
stant, .i.e., its energy density is almost constant and its 
equation of state (EOS) is close to w ~ — 1. On the 
other side, freezing fields have an EOS far away from the 
cosmological constant behavior, but that evolves towards 
w ~ — 1 at late times. 

More recently, there have been serious attempts to 
show that the dynamics of cosmological scalar fields 
should proceed under the action of some simple guidelines 
determined both by the scalar field equations of motion 
and the presence of other background fields needed for a 
complete description of the evolution of the Universe [5- 
[Toj . In a first approximation, we may think that the 
quintessence field is in a slow-roll regime similar to that 
of inflationary models, as in the two cases the main 
purpose is to provide of an accelerating expansion of 
the Universe. However, contrary to the inflationary 
case, the quintessence field, if really existent, has not 
yet dominated the expansion of the Universe and must 
have co-existed with other dominant components in the 
past. This co-existence should have made its imprint 
on the early evolution of the quintessence field and de- 
termined the present conditions for an accelerating Uni- 
verse. These simple arguments suffice to show that the 
quintessence dynamics cannot be simply described by the 
slow-roll formalism of inflation [8(. In fact, it has been 
shown that thawing and freezing fields have quite dis- 
tinct behaviors at early times. 



' lurcna@fisica.ugto.mx 



In this paper, we shall show that general guidelines ex- 
ist for the evolution of cosmological scalar fields, which 
are given in terms of the phase space structure of appro- 
priately chosen scalar field variables. These scalar field 
variables are simply the kinetic and potential contribu- 
tions of the scalar field to the total energy density of 
the Universe which, to our knowledge, were first defined 
in the seminal paper of Copeland, Liddle and Wands in 
Ref.[ll[ for the simple model of a quintessence field en- 
dowed with an exponential potential. 

Our claim is that dynamical attractor trajectories ex- 
ists in the form of the heteroclinic trajectories of the 
phase space, and that these special trajectories give an 
unified description and guidance for the dynamics of gen- 
eral scalar field models (see[I2,[l3| for a similar result for 
inflationary fields). We are then taking a step further 
from the standard approach in cosmological dynamical 
systems: it is not enough to search for the critical (fixed) 
points of the system, but also important is to determine 
the special trajectories that connect physically relevant 
critical points. In passing by, we shall argue that the 
phase space of the kinetic and potential variables is the 
appropriate arena for the description of the dynamics of 
general scalar field models, irrespective of the number of 
dimensions the true phase space needs to fully describe 
the evolution of the quintessence field. 

A brief description of the paper is as follows. In Sec.lIIl 
we present the general equations of motion and review 
the properties of the dynamical system corresponding to 
the case of an exponential scalar field potential; we will 
emphasize the role played by critical points and hetero- 
clinic lines in the determination of the scalar field dy- 
namics and the structure of the phase space. In Sec. IIIII 
we extend the description of the exponential potential to 
more general cases, again in terms of critical points and 
heteroclinic lines. In particular, we give a separate de- 
scription for the dynamics of thawing and freezing fields, 
and provide of simple explanations for the expected be- 
havior of quintessence fields. We present numerical re- 
sults in Sec. |IV] to support the qualitative description 
given in previous sections. Finally, Sec. [V] is devoted to 
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conclusions. 



and finally Eq. (ITal) reads 



II. DYNAMICAL SYSTEM 

Let us consider a flat Friedmann-Robertson- Walker 
(FRW) Universe that contains a perfect fluid with a 
barotropic equation of state (EOS) p 7 = (7— l)p 7 , where 
p 7 and p 7 are its pressure and its energy density, re- 
spectively. The known values of the EOS are 7 = 1 for 
pressureless matter (dust), and 7 = 4/3 for a relativistic 
fluid (radiation), but in general < 7 < 2. There is also 
a scalar field <j> which is endowed with a scalar potential 
V(4>). The calculations that follow will not depend upon 
the particular form V((f>), but we will focus our atten- 
tion on potentials that can show a useful quintessence 
behavior. 

The equations of motion, which together form the well 
known Einstein-Klein-Gordon (EKG) system of equa- 
tions, read 



H 



K 2 



p 7 = -3Hjp 7 , 



(la) 
(lb) 
(1c) 



Where H = a/ a is the Hubble parameter and a dot 
means derivative with respect to cosmic time. In ad- 
dition, one has the Friedmann constraint 
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Pi + ^ + v(4>) 



(2) 



with k 2 — 8ttG, and G is Newton's constant. The energy 
density of a homogeneous scalar field is p$ = {l/2)4> 2 + 
V{<f>), whereas its pressure is p^ = (l/2)0 2 — V{(f>). 

We follow here the dynamical system approach of 
Copeland, Liddle and Wands in Ref . fllj . If we define 
the new variables 



K(j> _ Ky/V 

~7M ' y= y/3H 



(3) 



then the Klein- Gordon (KG) equation of motion of the 
scalar field (jlc|) can be written as a dynamical system of 
the form 

x' = -3x + X^y 2 + ^x [2x 2 + 7 (1 - x 2 - y 2 )] (4a) 

V 1 = -^\[\ X V + \v i 2 * 2 + 7(1 - x 2 - y 2 )} , (4b) 

where we have defined the roll parameter A = 
—8^V/{kV), and a prime denotes a derivative with re- 
spect to the e- folding number TV = ln(a). The Friedmann 
constraint (|2|) now reads 



111 +x 2 +v 2 = 1 



(5) 



H 3 9 

B* = 2 ^ 



•7(1 



y 2 )] 



(6) 



Notice that it is not necessary to consider the equation 
of motion of the perfect fluid (jlbp . as we can use the 
Friedmann constraint ([5]) to determine p 7 once we have 
found a solution for x and y from Eqs. (j4}. We can 
also write the (non-constant) barotropic EOS 7^, and 
the density parameter f2^ of the scalar field, as 
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+ , 



2x 2 



Pd, 



3H 2 



x 2 + y 2 . (7) 



It is well known that the Eqs. (0} form an autonomous 
system if A = const., which is the case of an exponen- 
tial scalar potential of the form V(<fi) — Vbe~ AK *. Ex- 
ponential potentials are well understood dynamical sys- 
tems, as infered from the many papers that have stud- 
ied its dynamical properties in different cosmological 
settings[ll[ll[ll, see also 0] for a modern and gen- 
eral review on scalar field models. 

For a general scalar potential, it is necessary to con- 
sider an equation of motion for A itself, which in general 
terms is given by@, [H, 03 



A' 



-V6xX 2 (T- 1) 



(8) 



where T = Vd^V/ (d<j,V) 2 is the so-called tracker 
parameter 0. For those cases in which Eq. ((SJ can be 
written in terms only of variables x, y, A, then Eqs. (U) 
and ([5]) together form a three dimensional autonomous 
system. This was indeed the case thoroughly studied in 
Ref. QJ] . 

It is not our intention to do the same here, but we shall 
rather take a more general approach that includes even 
those cases in which Eq. ([SJ cannot be written in a closed 
form, and the equations of motion of higher derivatives 
of the potential would have to be taken into account. 



A. Critical points and stability 

We start by recalling the properties of the critical 
points, see also[l8j, (xo,yo)j those for which x'(xo,yo) — 
= y'(x 0} y ), of the 2-dim system ^ in the case of an 
exponential potential, A = const., which are shown in Ta- 
ble HI Critical points are also called fixed points, because 
they represent in themselves a solution of the dynami- 
cal system that appears on the phase space (x, y) as the 
single point (xo,yo). There are in general four types of 
critical points, whose existence and stability conditions 
depend upon the values of the perfect fluid EOS 7 and 
the roll parameter A. 

To study the stability properties of the critical points, 
one considers small perturbations of the form x = xq + u 
and y = yo + v, where (u,v) are the perturbation vari- 
ables. By taking a first order expansion of the full equa- 
tions of motion around a given critical point (xq, yo), 
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TABLE I. Critical points (xo, yo) of the dynamical system Q in the case of an exponential potential of the form V((j>) = Voe 
taken from Ref.[ll|. See text below for more details. 
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Stability 


"0 


7<p 


A 


n 

u 




n 


W A unn /*v 
V s\ dllU J 


kJdLlU.lt. JJU111L 1U1 U \ V \ 


n 


T In ri /"> 11 ri 

LJ llU.tllIlL.Ll 
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VA and 7 


Unstable node for A < \/E 
Saddle point for A > v6 


1 


2 


C 


-1 







VA and 7 


Unstable node for A > —y/6 
Saddle point for A < — VrJ 


l 


2 


D 






- A 2 /6 


A 2 < 6 


Stable node for A 2 < 37 
Saddle point for 37 < A 2 < 6 


1 


A 2 /3 


E 


^/372 7 /A 


^3(2- 


- 7)7/2A 2 


A 2 >3 7 


Stable node for 37 < A 2 < 24 7 2 /(97 - 2) 


3 7 /A 2 


7 



Stable spiral for A 2 > 24 7 2 /(97 - 2) 



the linear equations of motion can be written as 

The solution of the perturbation variables is formally 
given by u = e MN Ui, where the exponential term is itself 
a matrix that can be written in terms of the eigenvalues 
m and eigenvectors p, of the stability matrix M. Actu- 
ally, the general solution for a given critical point is of 
the form 

( " ) = Ci/2ie miW + C 2 fi 2 e m * N (10) 

where C\ , C 2 are arbitrary constants related to the initial 
conditions. 

If all eigenvalues have negative (positive) real part, 
Re(m) < 0(> 0), then the corresponding critical point is 
called stable (unstable); any imaginary part in the eigen- 
values will add an oscillatory feature to the solutions and 
they are further label as spirals. In the case there is one 
eigenvalue with a positive real part and another one with 
a negative real part, then the critical point is called a 
saddle. 

For each one of the eigenvalues there is one eigenvec- 
tor, and hence there are two eigenvectors for each critical 
point. Each eigenvector jl determines a principal direc- 
tion around a critical point along which the solution is 
purely given by the corresponding eigenfunction e mN . 



B. Phase space structure 

Before we proceed further, we analyze the structure of 
the phase space in the case of an exponential potential. It 
is usually assumed that one should only take care of the 
upper part of phase space y > for expanding Universes. 
However, the fact is that the upper and lower parts of the 
phase space are separated regions, each one enclosed by 
boder lines called as separatrix. In more strict terms, a 
separatrix is not but the heteroclinic line that connnects 
two different critical points. 



Heteroclinic lines depart from an unstable or saddle 
point along its unstable principal direction, and arrives 
at another stable or saddle point along its stable principal 
direction. In the case of Eqs. (U), we can identify four 
or five heteroclinic lines, depending on the values of A, 
as the latter determines the number and nature of the 
critical points, see Table HI 

One can show that there are three heteroclinic tra- 
jectories that correspond to the conditions y — and 
x 2 + y 2 = 1 (it can be easily verified by direct substi- 
tution that these conditions are exact solutions of the 
full dynamical system Q). The heteroclinic trajectories 
for (x, y = 0) are those that depart from the unstable 
kinetic-dominated points (±1,0) (B and C, respectively, 
in Table HJ) and arrive at the saddle perfect fluid domi- 
nated point (0, 0) (point A in Table[I]). These heteroclinic 
lines are labeled as II and III in the plots in Fig.[TJ and are 
also the separatrices of the y > and y < regions of the 
2-dim phase space. In the case x 2 + y 2 = 1 , we can have 
one heteroclinic trajectory departing from the unstable 
kinetic point (—1,0) and arriving to the the saddle kinetic 
point (1,0) if A 2 > 3 7 (line IV in bottom plot in Fig. [J); 
or we can have two heteroclinic lines departing each one 
from any of the unstable kinetic points (±1,0) and arriv- 
ing to the (saddle or stable) scalar field dominated point, 
corresponding to 3 7 < A 2 < 6 or A 2 < 3 7 , respectively 
(lines IV and V in top and middle plots in Fig. [1]). In 
any case, the heteroclinic lines x 2 + y 2 — 1 do not allow 
variables x, y to take values larger than unity. This is 
a manifestation of spatial flatness as imposed upon our 
model by the Friedmann constraint ([5])- 

There is not an analytic expression for the hetero- 
clinic line that connects the saddle perfect fluid dom- 
inated point (0,0) (A in Table H]) to the stable scalar 
field-dominated solution (for A 2 < 3 7 , point D in Ta- 
ble H| or to the scaling solution (for A 2 > 3 7 , point A in 
Table IJ), that we have labeled as I in all cases in Fig. [1] 
There is also one more heteroclinic line that joins the 
saddle scalar field-dominated and the stable scaling so- 
lutions, but it only appears for 3 7 < A 2 < 6 (line VI in 
middle plot in Fig. [1]) . 

As we said before, heteroclinic lines split the phase 
space in separated regions that are disconnected one from 
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each other, that is, any trajectory must remain within 
the region it initially started in. There are two separated 
regions for A 2 < 6, but only one single bounded region if 
A 2 > 6. 



C. General dynamics for the exponential case 

The general dynamics of a scalar field endowed with 
an exponential potential is completely described by the 
properties of the critical points and the heteroclinic lines 



of its phase space. Examples of different trajectories were 
given m Rcf.lllL Here we will focus our attention only in 
trajectories that are consistent with an early perfect-fluid 
dominated Universe. 

In Fig. [T]we also show examples of different trajecto- 
ries departing closely from the perfect fluid domination 
point A. The general behavior of these particular trajec- 
tories can be quoted quite simply: each trajectory lies 
close to the heteroclinic line that joins the saddle perfect 
fluid-dominated point A and the stable critical point at 
hand in each case (see the heteroclinic lines I in all plots 
in Fig. [1]) . We can turn this quote into a stronger state- 
ment: the general dynamics of a cosmological scalar field 
endowed with an exponential potential, is well described 
by the heteroclinic trajectory that joins the perfect fluid 
dominated point with the late time attractor present in 
the phase space. 

By studying this single trajectory, one can in princi- 
ple determine the agreement of a given model, in the 
exponential case, with the expected evolution of the Uni- 
verse. From the numerical point of view, the drawing of 
the heteroclinic trajectory needs of very especial initial 
conditions. It is not possible to find a full analytical so- 
lution for the heteroclinic trajectory we are interested in, 
but at least the behavior of the trajectory close to the 
origin of coordinates can be easily determined as follows. 

We first write the linear equations of motion that 
are used to study the stability of critical points in au- 
tonomous systems. For the case of the perfect fluid- 
dominated point A, the linearization of the system (U) 
has the explicit solution 



e -3(2- 7 )jV/2 + 



Vi 



(11) 



where we used the original phase space variables (x, y) 
for the perturbations (u,v), which is very convenient as 
the perturbation procedure is made around the origin of 
coordinates. Here, Xi and yi refer to the initial values of 
the phase space variables. 

For both cases of a relativistic and a pressureless fluid, 
the first eigenvalue is negative definite, m\ = —3(2 — 
7)/2 < 0, and then the kinetic energy of the scalar field 
must decrease for trajectories nearby the critical point 
A. On the other hand, the second eigenvalue is positive 
definite, m% = 37/2 > 0, and then the potential energy 
must increase. In terms of the dynamical system jargon, 
one says that the critical point A is a saddle point and 
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FIG. 1. Structure of the phase space according to the dynam- 
ical system Q in the case of A = const. We show the critical 
points of Tableland the heteroclinic trajectories described in 
Sec.[IT]that j°i n them in each case, see text for details. Notice 
that the phase space is split in different disjoint regions for 
which the heteroclinic lines act as separatrices. Also shown 
are different trajectories in phase space obtained as numeri- 
cal solutions of the dynamical system @, that depart from 
different points with yt = 10 -4 . It can be seen that all trajec- 
tories are focused around the heteroclinic line I that connects 
the saddle point A (see the insets) with the corresponding 
stable point D or E in each case. The semicircle around the 
origin of coordinates corresponds to Q.^ = x 2 + y 2 = (0.05) 2 . 
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its normal directions, represented by the eigenvectors of 
the stability matrix A4, see Eq. ([3]), coincide with those 
of the phase space axes, see Eq. (fTl"]) . 

We have learned then that the kinetic energy of the 
scalar field decreases if the latter is subdominant with 
respect to the perfect fluid, which is the expected situ- 
ation in the early Universe. This also means that the 
kinetic variable x must be smaller than the potential one 
y. If we now make a second order expansion of Eq. (I4al) . 
but still a first order expansion of Eq. (lib)) , around the 
same critical point A, we find 



3(2-7) , . /3 2 
/ 3 7 



y 



-y . 



(12a) 
(12b) 



whose general solutions, under the assumption that A is 
constant, are 



y(N)= yi e^ N / 2 . 



2 A 



3 (2 + 7 ) 



V 



(13a) 
(13b) 



It can be seen that, if we neglect the decaying solution, 
the kinetic energy evolves proportionally to the square of 
the potential energy, namely 



x{y) 



X 



(2 + 7) V6 



y 



(14) 



This is in good agreement with our assumption that the 
kinetic perturbation x is one order of magnitude smaller 
than potential perturbation y. For future reference, we 
shall call Eq. (|14|) the thawing constraint. Notice that 
Eq. (fT4| has the same form as the famous slow-roll con- 
dition 



3H(/) ~ -dsV 



x{y) 



A 

7! 



y 



(15) 



except for the correction induced by the presence of the 
perfect fluid through its EOS 7. 

We have found the early behavior of the heteroclinic 
line that departs from the perfect fluid-dominated point, 
in the phase space the trajectory corresponds to the loci 
of a rotated parabola. Such a behavior can be easily seen 
in all lines I in Fig. [TJ actually, the heteroclinic trajecto- 
ries in Fig. [T] were found from the numerical solution of 
Eqs. (U]) with initial conditions subjected to Eq. (fT4"|). 

Eq. (|14[) can be written in the more recognizable terms 
of the flow parameter F = r y ( p/(fl ( f ) X 2 ), defined by Cahn, 
de Putter, and Linder in Rcf.|8]. If we take the thawing 
constraint (fT4"]) and the definitions in Eq. (J7J, we find 
that 



F 



3(2 + 7 ) 2 



(16) 



and then we recover the known value F = 4/27 if the 
dominant perfect fluid is pressureless matter (7 = 1), 



whereas F = 3/25 if it is a relativistic fluid (7 = 4/3). In 
passing by, we note that the value of the flow parameter 
under the terms of the slow- roll condition (fT5|) is F = 1/3. 
The asymptotic values of some other quantities around 
the point (0,0) can also be found, as for instance d, [1]: 



3H<j) 
~d A V 



d<t,V 



2 + 7 



(17) 



As we have said before, critical points in the exponen- 
tial case are also fixed points, and then they seem to have 
more relevance that the trajectories themselves on phase 
space. This is true in the case of early and late time 
attractors, but particular trajectories are important if in 
addition one is interested in solutions which are in agree- 
ment with the stringent constraints arising from the ra- 
diation and matter dominated eras of the early Universe. 

In Fig. [1] we show that trajectories with initial val- 
ues close to the perfect fluid-dominated point A follow a 
quite similar path to that of the heteroclinic trajectory 
that connects point A with any of the stable fixed points 
at hand. The seemingly attractor character, or focusing 
property, of this heteroclinic line comes from a combina- 
tion of the saddle properties of point A and the attractor 
character of the ending stable point. In the following sec- 
tions, we will argue that this nice property is preserved 
in the case of general scalar field potentials. 



III. DYNAMICS OF GENERAL SCALAR FIELD 
MODELS 



In the case of an exponential potential, critical points 
are truly equilibrium fixed points. For the case A is also 
a variable, those points are not fixed, but they still repre- 
sent, for a given value of A, the points at which the 2-dim 
phase space velocity (x' , y') vanishes. In other words, the 
velocity field on the 2-dim phase space (x, y) has the same 
structure than that in the exponential case, but now the 
location of the critical points can change as A evolves. 

If A changes slowly, the critical points can be treated 
as if they were approximately fixed points; this feature 
is particularly useful in the case of quintessence fields, 
as one expects slow variations in the roll parameter A 
to have a solution with an accelerating expansion of the 
Universe at late times when the scalar field is the domi- 
nant component. However, the same can happen at early 
times, when the scalar field is subdominant with respect 
to the radiation and matter components. Hence, we ex- 
pect the early and late time evolution of general scalar 
field models to proceed as in the case of the exponential 
potential, even though the intermediate-time dynamics 
of the scalar field may be more involved due to the evo- 
lution of the roll parameter A. 
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A. Early dynamics of thawing and freezing scalar 
field solutions 

From the dynamical system perspective, A plays an 
important role in determining the existence and stability 
properties of two important critical points: the scalar 
field dominated solution D and the scaling solution E, 
see Table HI As we shall show, some aspects of the early 
and late time evolution of scalar fields can be understood 
if those two critical points are taken into account. 

For simplicity, let us start with the assumption that 
A -C 37, and that A is a monotonic growing function of 
N (this condition is not crucial for the argument, but it 
helps to visualize the description below). Thus, the scal- 
ing solution does not exist, and the only critical solution 
on the 2-dim phase space is that of scalar field domina- 
tion, point D. If A changes slowly and keeps a value be- 
low 37, the 2-dim velocity flow will focus the trajectories 
to that point, and then we shall recover a scalar field- 
dominated and accelerating Universe as long as A < v2 
at the present time. 

Notice that wc have just described the evolution ex- 
pected for thawing solutions, in which the EOS starts 
very close to the cosmological constant value, 7^ ~ 0, 
the flow parameter is F = 4/27, and the field gains ki- 
netic energy as its evolution proceeds in the form sug- 
gested by Eqs. (TT4")) . As the value of A grows, the scalar 
field-dominated point D moves clockwise and so does the 
phase space trajectory (see, for instance, Fig. [5] below). 

Another extreme case is that in which A 2 3> 37, and 
A is a monotonic decreasing function of N (as before, 
this condition is not crucial for the argument, but it also 
helps to visualize the description below). For this case, 
the scalar field-dominated solution D does not exist, but 
the scaling solution E comes into play. For large val- 
ues of A, this latter critical point is located very close to 
the perfect fluid domination one, point A, see Table fl] 
Again, if the variation of A is slow enough, the field fol- 
lows the heteroclinic trajectory that departs from the 
perfect fluid-dominated solution, hooks up with the scal- 
ing solution and follows it as the value of A changes. It 
should be stressed out here that we are referring to the 
same heteroclinic line described above for thawing fields 
and which is partially described by Eq. (fT4")) . 

The scaling point E moves outwards from the origin in 
phase space if A decreases, and then 7^ ~ 7 throughout 
this period; if at some point A 2 < 37, then the scaling 
solution disappears and the scalar field then looks for 
the scalar field-dominated solution D. Contrary to the 
previous case, the trajectory of the scalar field moves 
counterclockwise in the phase space (see, for instance, 
Fig. [3] below) . While the scalar field sits down on the 
scaling solution, we easily infer that the flow parameter 
is F — 1/3, and then we have just described the evolution 
of freezing tracker fields. 

We have seen that for a tracker solution we need the 
value of A to allow the appearance of the scaling solution 
close to the perfect fluid dominating solution, which is 



assured by asking A 2 3> 37. For intermediate values A ~ 
37, the scaling solution may not exist or, if it exists, it 
is located far away from the origin of the phase space. 
Thus, the behavior in the proximity of the origin must 
be that of the thawing solution. However, notice that a 
purely thawing potential, for which A grows, will never 
be able to accelerate the expansion of the Universe at 
late times if initially A ~ 37. 

If the same is assumed for a purely freezing potential, 
for which A decreases, then the scalar field will anyway 
show a thawing behavior, because the scaling solution 
is never at hand for it to catch it up. The freezing po- 
tential has the chance to accelerate the expansion of the 
Universe, if A decreases fast enough as to fulfill the con- 
dition A < \/2 before the present time. 



B. Geometric interpretation of the Flow parameter 

There is a simple geometrical interpretation of the flow 
parameter F in Eq. (|16|) in terms of trajectories in phase 
space. As we have said before, what we have called the 
thawing solution (| 14[) indeed corresponds to the initial 
part of the heteroclinic trajectory that joins the (saddle) 
perfect fluid-dominated solution A to any of the stable 
solutions at hand. Most trajectories seem then to be 
focused towards the heteroclinic trajectory, and then the 
adjective flow is well deserved for the parameter F. 

In the case of freezing behavior, we have shown that 
trajectories first engage with the thawing solution and 
then hook up with the scaling solution; in consequence, 
the flow parameter takes on the freezing value F = 1/3. 
The flow now refers to the fact that solutions are follow- 
ing the scaling solution as it evolves outwards from the 
phase space origin. The focusing of the trajectories is 
achieved via the stability properties of the scaling solu- 
tion. 



C. The tracker theorem 

The nice attractor features of freezing fields at early 
times come basically from its tracker nature, a term first 
applied to scalar fields in Ref.Q. In this paper, the 
authors established a theorem for the necessary condi- 
tions for a scalar field to show tracking behavior; in sim- 
ple terms, one needs at early times the tracker param- 
eter to be nearly constant and to satisfy the condition 
(r — 1) > 0. One can though understand the tracker 
conditions in more simple terms by using the dynamical 
system equations ((4]). 

We have seen that the attractor property of freezing 
fields arises because one sets up large initial values of 
the roll parameter A, so that the scaling solution E is a 
critical point of the 2-dim phase space velocity (x',y'). 
The larger the value of A, the easier is for the scalar 
field to catch up the scaling solution. From this point 
of view, the tracker condition just refers to large initial 
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values of the roll parameter that allow the appearance of 
the scaling solution E. (As a curiosity, one can see that 
the attractor character of the tracker solution discussed 
in Ref.0] actually coincides with the analysis made in 
Ref.pH for the scaling solution.) 

However, one can go back to the usual tracker condi- 
tions if we take a look at the evolution equation of the 
roll parameter, Eq. <JSj> - It can be seen that large initial 
values of the roll parameter are only allowed for scalar 
field models for which A is a growing function if we go 
backwards in time; that is, we need (r— 1) > 0. What we 
do not need is the slow variation of the tracker parameter, 
because, as we shall show in the numerical experiments of 
Sec.|lVl the very condition A> 1 leads to x — > 0. Hence, 
it is not the tracker parameter but the vanishing of the 
kinetic energy the responsible for the slow-variation of 
the roll parameter at early times, see Eq. ((HJ. This en- 
hances the possibilities of the scalar field to catch up the 
scaling solution, as the corresponding critical point on 
the 2-dim phase space (x, y) moves adiabatically. 

Therefore, the tracker theorem can be restated in sim- 
pler terms as follows. The tracker property appears for 
any scalar field model in which the roll parameter A is 
capable of taking on large initial values in the early Uni- 
verse. It is easier to look at the expected evolution of A, 
than to explore the more involved tracker parameter T, 
to determine the tracker properties of a given scalar field 
model. 



D. Late time dynamics 

The late time dynamics of the quintessence fields will 
depend upon its thawing or freezing behavior. In the 
case of thawing fields, we expect that they slowly deviate 
from the cosmological constant values, and that A moves 
from small to large values. If a thawing field is currently 
accelerating the expansion of the Universe, then A < \J~2 
today. It may be the case that A takes even larger values 
in the future, so that the quintessence field is not able 
to accelerate the cosmic expansion and can even provoke 
the appearance again of the scaling solution E on the 
phase space if at some point A 2 > 37. 

As for freezing fields, they start with large values of A, 
but the latter has to migrate to smaller values A < \/2 
so that the inflationary scalar field-dominated solution 
D appears on time. If this were also the trend for the 
future, we can only expect that a freezing field will only 
come closer to the cosmological constant behavior, i.e. 
x — s- and y —> 1, as A — > 0. 

Whether the field is thawing or freezing, its capability 
for providing an accelerating solution nowadays can be 
characterized in terms of A, as it is only necessary that the 
(inflationary) scalar field dominated solution D becomes 
alive. Hence, nowadays A should be of the order of A 2 ~ 
3(1 + w^o). 

Even if the early dynamics of a scalar field resembles 
very much that of an exponential potential, the late-time 



dynamics can be more involved and strongly depends 
upon the explicit form of the scalar field potential (Bl fl6|. 
Nonetheless, an easy description can be given for those 
cases in which Eqs. (UJ) and ((SJ together form a closed and 
autonomous system of equations [16|. For these cases, the 
phase space is truly 3-dimensional, (x, y, A), and we will 
refer to it hereafter as the augmented phase space. No- 
tice that the heteroclinic lines I, II, and III described 
in Sec. Ill Bl retain its character for any value of A; in 
geometric terms, we can see that the augmented phase 
space occupies, in general, the interior region of a unitary 
cylinder whose axis lies along the A-direction. 

If there is a value A* such that [r(A*) — 1] = 0, then 
there exist generalized scalar field-dominated and scaling 
solutions in the augmented phase space (x, y, A). The 
properties of these critical fixed points are quite similar 
to their exponential counterparts in Table U in that one 
only needs to make the replacement A — > A*. Moreover, 
there also exists the generalized perfect fluid-dominated 
solution, which exists for any value A and is always a 
saddle point. 

Whenever the generalized scalar field dominated D and 
scaling E solutions exist, there can be a unified descrip- 
tion of the dynamics in terms again of the heteroclinic 
line that joins the perfect fluid dominated point A to the 
stable point available in the augmented phase space. As 
in the case of the exponential potential, the general solu- 
tion for a given model is then characterized by such het- 
eroclinic line. Even more, in our numerical experiments 
with different potentials in Sec. IIV1 we will find that the 
early description of the heteroclinic points is also given 
by the thawing behavior in Eq. (fiH) . This was to be ex- 
pected, as the evolution of A depends upon the second 
order perturbations of x, see Eq. ([8]). and then most of 
the motion around the perfect fluid-dominated solution 
must proceed at A ~ const. 

There is though a non-trivial difference with respect to 
the exponential case: the de Sitter point (0, 1) exists as 
a late time attractor if A — > and [r(0) — 1] > 0. That 
is, for any scalar field model in which the roll parameter 
vanishes at late times the final fate is the cosmological 
constant easeful Il6l]. 



IV. NUMERICAL EXAMPLES 

We now present some numerical solutions of the EKG 
equations in cases in which Eqs. Q and ((8j) form a closed 
and autonomous system of equations. This is done for 
numerical purposes, as in other more general cases we 
may need to write a full hierarchy of equations for higher 
derivatives of the scalar field potential. As we have dis- 
cussed before, it will be only necessary to look at the 
properties of the 2-dim phase space (x, y) to have an am- 
ple picture about the dynamics of scalar fields. 
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A,=100 




FIG. 2. Numerical solutions of Eqs. @ and ((IB) , for nu- 
merical convenience we chose n = 50. Shown are the tra- 
jectories (black-solid and blue-dotted lines) in the 3D phase 
space (x, y, A) for different initial conditions of the system ((J) 
and ()18|) . The corresponding projections onto the (x,y) plane 
are also shown (red-solid and green-dashed lines); see the text 
and Fig. [3] for a comprehensive description of the dynamics 
on the rry-plane. 



A. Monotonic freezing case 

As a typical example for a freezing model we take the 
tracker models V(<j») = M 4+n cj)~ n , with n > 0, for which 
the roll parameter and the tracker parameter are, respec- 
tively, A = n/ (k<P) and T— l = l/n>0; notice that this 
is a monotonic freezing case, as the latter condition on 
the tracker parameter means that the roll parameter is 
an always decreasing function. In fact, the equation of 
motion of the roll parameter is 



A' 



V6 2 

x\ . 



(18) 



and then it is one of those fortunate cases in which the 
equations of motion Q and (|18[) together form an au- 
tonomous system of equations, which is able to represent 
the full dynamics of the original EKG system. 

According to the phase space analysis in previous sec- 
tions, this type of models have the following fixed points 
in the augmented phase space on the plane A = 0: per- 
fect fluid-dominated point A, kinetic dominated points B 
and C, and the de Sitter point (0, 1). In addition, points 
D and E are also critical, but not fixed, points at differ- 
ent slices with A = const. Even th oug h we cannot apply 
here the stability analysis of Ref . |16fl . we expect the de 
Sitter point to be a stable point: once the field is rolling 
down its potential we have x > 0, and then A should 
be an ever decreasing function according to Eq. (|T8|) . In 
this form, the system must approach the de Sitter point 
asymptotically. 

Numerical solutions of the set of Eqs. (0} and (|18p are 
given in Figs. [2] for the 3-dim phase space (x, y, A), and 
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FIG. 3. Numerical solutions of Eqs. © and (fig)) , which to- 
gether form an autonomous dynamical system; for numerical 
convenience we chose n = 50. (Top) Trajectories of the nu- 
merical solutions on the 2-dim (x,y) plane. The dots denote 
the positions of the critical points D and E, see Table [I] as 
A evolves. All the solutions initially have Ai = 100. In par- 
ticular, the (red solid) line is the heteroclinic line that joins 
the perfect fluid-dominated point A with the de Sitter point 
(0, 1). (Bottom) The evolution of the scalar field EOS w$ 
(red-solid and green-dashed lines) and the flow parameter F 
(black-solid and blue-dotted lines) for the same numerical so- 
lutions as in the top figure. Both parameters show a tracker 
intermediate behavior for which — const, and F = 1/3. 
At late time, the EOS takes on negative values as the solution 
evolves towards the de Sitter point and enters a pure slow-roll 
regime with F — > 1/3 again, see Eq. (|15|) . 



in Fig. [3] for the 2-dim phase space (x, y), for n = 50. A 
detailed description of the behavior of the different curves 
is as follows. The attractor trajectory is the heteroclinic 
line that joins the saddle point A (0, 0) with the stable de 
Sitter point (0, 1 the behavior of this heteroclinic line 
in the neighborhood of point A is given by the thawing 



1 We refer here to a 2-dim heteroclinic line that connects the criti- 
cal points (0, 0) and (0, 1), even though the true heteroclinic line 
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solution (fT4|) . Because of the large initial value of the 
roll parameter, A^ = 100, the scaling solution, point E, 
is present. The heteroclinic trajectory hooks up with the 
scaling point, and the latter moves outwards as the value 
of A decreases. Once A < V6, there also appears the 
scalar field dominated solution, point D, which starts to 
move counterclockwise on the unitary circumference as 
A continues decreasing. Both critical points merge once 
A = 37 and the scaling point formally disappears, and the 
motion of the heteroclinic line is now guided by the still 
moving point D until they reach together the de Sitter 
point. 

The description we have given for the heteroclinic line 
suffices to understand the general behavior of inverse 
power-law potentials. To show the attractor feature of 
the heteroclinic line, we have also solved the equations 
of motion for diverse initial conditions (xi,yi) but keep- 
ing the same A^ = 100. It is clear that the curves follow 
the heteroclinic line at late times, as expected. But for 
this type of potentials we notice a peculiar early behav- 
ior denoted by the clockwise semi-circumferences. This 
behavior is easily explained if we write the equations of 
motion (Q} in the limit A ^> 1, and then 




It can be shown that, under these equations, the scalar 
field density parameter ([7|) is a conserved quantity, i.e. 
CI,), — x 2 + y 2 — const., and then the phase space vari- 
ables should move clockwise on circumferences until they 
lose enough potential energy (y — > 0) and the scalar 
field energy is dominated by its kinetic component (there 
is no other option, anyway, as the variable y cannot 
become negative because of the heteroclinic separatrix 
(x,y = 0)). The clockwise motion is explained by the 
fact that the field prefers to move downhill its potential, 
x > 0, rather than otherwise. We finally note that the 
circular motion occurs quite rapidly, an estimation can 
be made from Eqs. (fH?|) . and then the e-fold interval for 
the circular motion is AiV ~ 1/A. 

In Fig. [3] we also make a translation of the scalar field 
dynamics in terms of the scalar field EOS = 70— 1 (see 
Eq. 0) and the flow parameter F (see Eq. ((TS])). The 
scalar field EOS reaches a tracker regime that almost per- 
fectly traces the background EOS, which is 7 = 1 for the 
numerical examples; this is actually the tracker behavior 
described in detail in Ref.[2]. The small difference be- 
tween the two EOS means that the scalar field redshifts 
a bit slower than the background field. This difference, 
even if small, is important for the later overtaking of the 
background dynamics by the scalar field. The value of w$ 
asymptotically approaches the de Sitter value, w^, — » — 1. 



exists in the augmented 3-dim space (x,y, A), and connects the 
points (0,0,0) and (0, 1,0)Q]J, see also Fig. E] We will use the 
2-dim nomenclature as it fits well what we want to describe on 
the 2-dim phase space. 



As for the flow parameter, its value in the tracker 
regime is the expected one, F = 1/3, which is a con- 
stant value for almost the whole evolution. There is 
small 'bump' at the transition point at which the solution 
leaves the tracker regime, but F^ — > 1/3 asymptotically 
again, signaling that the scalar field is entering into a 
pure slow- roll regime at late times, see Eq. (fT5j) . 

As the heteroclinic trajectory is a strong attractor, it 
suffices to scan the form of heteroclinic trajectories in 
phase space for different initial values of A in order to find 
physically acceptable accelerating solutions, as is made in 
Fig. [3J We can see that the initial value of A should be 
small enough so that the heteroclinic line can provide 
of physical parameters in agreement with cosmological 
observations. According to Fig. HI it is necessary to have 
Xi < 0.5 if we want heteroclinic trajectories to be inside 
the arc of radius £1$ = 0.7 and angles corresponding to 
-1 < u> < -0.94. 
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FIG. 4. Heteroclinic lines departing from the perfect fluid- 
dominated point A, that were obtained from the numerical 
solutions of Eqs. ((J) and ()18[) . For each case, the initial con- 
ditions for x and y were chosen according to the thawing 
constraint (|13p . The shaded region is a circular arc of radius 
Q,^ = 0.7 and angles corresponding to — 1 < to^ < —0.94. It 
can be seen that only heteroclinic lines with Ai < 0.5 may be 
able to provide of a physically viable accelerating Universe. 



B. Monotonic thawing case 

To study the general behavior of thawing models, we 
may use the example of a PNGB potential of the form 
V{4>) = M 4 cos 2 (amp) , where a is a free parameter of the 
modelQ. The roll and tracker parameters for this kind 
of potentials are 



A = 2atan(aK0) , T(A) = — 




The PNGB potential is a representative case of a mono- 
tonic thawing case, for which the roll parameter is al- 
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FIG. 5. Numerical solutions of Eqs. @ and (JTSJ, for nu- 
merical convenience we chose n — 20. Shown are the trajec- 
tories (black-solid and blue-dotted lines) in the 3-dim phase 
space (x, y, A) for different initial conditions of the system @ 
and (|21|l . The corresponding projections onto the (x, y) plane 
are also shown (red-solid and green-dashed lines); see the text 
and Fig. [6] for a comprehensive description of the dynamics 
on the xy-plane. 



ways an increasing function. However, we should re- 
member that the PNGB potential has a zero minimum 
at aK(f> — 7r/2, and it is at this point that the roll param- 
eter diverges, A — > oo; hence, the autonomous dynamical 
system cannot describe the dynamics of the scalar field 
model beyond this point. This is a caveat of the dynam- 
ical system approach that appears for any model with a 
zero minimum, see for instance (l3l. Il9j . 

Some other models without a minimum in the poten- 
tial can be chosen from the examples in the specialized 
literature. We shall use a generic example in which the 
representative features of monotonic thawing potentials 
appear acutely. We write Eq. (JSJ in the form 



V \/6 2 
A = x\ 



(21) 



which also arises in the case of power law potentials 
V ~ 4> n . We have verified that other monotonic thawing 
potentials display qualitatively the same generic behavior 
of this example (see@ for some other instances). 

As in the case of freezing models, thawing potentials 
has the following fixed points in the augmented phase 
space on the plane A = 0: perfect fluid-dominated point 
A, kinetic dominated points B and C, and the de Sitter 
point (0, 1). In addition, points D and E are also critical, 
but not fixed, points at different slices with A = const. 
Contrary to the freezing case, we expect the de Sitter 
point to be an unstable point: once the field is rolling 
down its potential we have x > 0, and then A should be 
an ever increasing function according to Eq. ([21]) . This 
latter fact has another undesirable consequence for thaw- 
ing models: there is not a fixed critical point that can 
act as a late-time attractor. Thus, the thawing solution 



0.5 | 



5° o >y 




FIG. 6. Numerical solutions of Eqs. Q and (|2T|) . which to- 
gether form an autonomous dynamical system; for numerical 
convenience we chose n = 20. (Top) Trajectories of the nu- 
merical solutions on the 2-dim (x,y) plane. The dots denote 
the positions of the critical points D and E, see Table [2 as A 
evolves. All the solutions initially have Ai = 2. In particular, 
the (red solid) line is the homoclinic line that joins the perfect 
fluid-dominated point A with itself. (Bottom) The evolution 
of the scalar field EOS (red-solid and green-dashed lines) 
and the flow parameter F (black-solid and blue-dotted lines) 
for the same numerical solutions as in the top figure. Both pa- 
rameters show an early thawing behavior for which w$ ~ — 1 
and F = 4/27. At late times, the EOS takes on positive values 
as the solution evolves back to the point (0, 0) and F — > 1/3. 



we have found before in the neighborhood of the perfect 
fluid-dominated solution, see Eq. (|13j) . is no longer part 
of a heteroclinic line, but now we have the presence of a 
homoclinic line only, that eventually connects a critical 
point with itself. 

Numerical solutions of the set of Eqs. (0} and (|21l) are 
given in Fig. [5] for the 3-dim phase space (x,y,X), and 
in Fig. |6] for the 2-dim phase space (x,y); in all cases 
Ai = 2. The same features appear in both Figs. 02 and0 
so that a quick comparison can be made of the thawing 
case with the freezing case. 

First of all, we cannot notice any kind of strong at- 
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tractor trajectory on the phase space at early times as 
it is the case in freezing models. This was anticipated 
above because there is only a homoclinic trajectory de- 
parting from and eventually arriving at the (saddle) per- 
fect fluid-dominated point; in practice, the homoclinic 
line is a closed loop0. The departing behavior of the loop 
is given by the thawing condition (|13[) . but the late time 
behavior is dictated by the scaling solution E as A — > oo. 
The homoclinic trajectory, however, do split the (x, y) in 
two parts, as trajectories departing at x < (x > 0) sur- 
round the loop in the clockwise direction by its outside 
(inside) part as they approach the scaling solution E at 
late times. 

In terms of physical variables, we notice that the scalar 
field EOS for the homoclinic line departs from — — 1, 
and ends up with a constant value at late times which is 
a bit stifier than the background EOS (in this example 
7 = 1). As for the flow parameter, we notice the expected 
thawing behavior at early times, F = 4/27, but also the 
freezing one at late times, F = 1/3. 

Even though the homoclinic trajectory is not a strong 
attractor, it anyway determines the way other arbitrary 
trajectories wander about on the phase space. Thus, in 
order to search for valid accelerating solutions, it suffices 
to scan the form of homoclinic trajectories in phase space 
for different initial values A^. This is actually shown in 
Fig. [7J For small enough values of A, the homoclinic 
trajectory reaches the unitary circle because of the pres- 
ence of the scalar field-dominated point D. Remember 
that this critical point exists for A 2 < 6 and is stable for 
A 2 < 2. We also show in Fig. [7] the homoclinic trajecto- 
ries that remain within the arc of radius f2^ = 0.7 and 
angles corresponding to —1 < < 0.94; these trajec- 
tories could provide a valid accelerating Universe at the 
present time. 

V. CONCLUSIONS 

We have given evidence for the existence of an uni- 
fied description for the dynamics of general scalar field 
models, in terms of the heteroclinic trajectory that con- 
nects the perfect fluid-dominated stage of the Universe 
with any of the late-time attractors at hand in the phase 
space of the scalar field variables. 

The numerical results shown in Figs. |4] and [7] strongly 
suggest that, in the case of single field models with a 
monotonic evolution of the roll parameter A, the field 
should start in a thawing regime far away from the 
tracker one if it is to be consistent with an accelerat- 
ing Universe at late times. Actually, the tracker regime 



2 In strict sense, what we have called a homoclinic trajectory on 
the (x, y) plane is not a loop, but actually an open line in the 
augmented 3-dim phase space (x, y, A). It departs from the point 
(0, 0, 0) and arrives to the point (0, 0, oo), see for instance Fig. 131 
We shall use the term homoclinic anyway as it encloses well what 
we see on the (x, y) plane 
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FIG. 7. Homoclinic lines departing from the perfect fluid- 
dominated point A, that were obtained from the numerical 
solutions of Eqs. Q and (|21[) . For each case, the initial con- 
ditions for x and y were chosen according to the thawing 
constraint (|13[) . The shaded region is a circular arc of radius 
= 0.7 and angles corresponding to — 1 < to^ < —0.94. It 
can be seen that only homoclinic lines with Xi < 0.5 may be 
able to provide of a physically viable accelerating Universe. 

seems to spoil in all cases the viability of quintessence 
models even if its early-time attractor properties are very 
appealing |2.l20j. 

However, the early thawing behavior is a problematic 
one, as it requires special initial conditions, in particu- 
lar, small initial values of the roll parameter A. More- 
over, as in the case of an exponential potential [ll|, the 
initial value of the quintessence parameter fi^ must be 
extremely small in order to be consistent with cosmolog- 
ical constraints and to prevent quintessence domination 
well before the present time. As for the late time dy- 
namics, the approach presented here suggests that the 
roll parameter A is small at late times, and that its value 
should be related to the present value of the dark energy 
EOS by A 2 ~ 3(1 + wq); in other words, the derivative of 
the quintessence potential must be related to the actual 
value of the dark energy EOS. 

It is interesting to note that the aforementioned diffi- 
culties for thawing and freezing scalar fields are act ually 
well known from the case of an exponential potential}! lj. 
They are actually present in more general cases because 
the phase space equations of motion preserve its form in 
the general case, and then quintessence fields inherit the 
difficulties of the exponential case. 

There is still the possibility of having single field mod- 
els with a non-monotonic evolution of the roll parameter, 
in which the strong attractor properties of freezing fields 
can be combined with the appropriate late-time dynam- 
ics of thawing fields. This can in principle be achieved 
with the use of very contrived potentials (see an early 
example in[2|), in the case in which the roll parame- 
ter changes stochastically, or even in the case in which 
phase space dynamics needs a full hierarchy of equations 
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of motion [l7}. 

In any case, our general suggestion in the study of 
scalar field models is still the same. One should look 
for the heteroclinic line departing from the perfect fluid- 
domination point; this single trajectory encloses in itself 
enough information to determine the capabilities of the 
given quintessence field to act as a realistic dark energy 
model. 
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